Second-Order Cosmological Perturbations from Inflation 
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We present the first computation of the cosmological perturbations generated during inflation 
up to second order in deviations from the homogeneous background solution. Our results, which 
fuUy account for the inflaton self-interactions as well as for the second-order fluctuations of the 
background metric, provide the exact expression for the gauge-invariant curvature perturbation 
bispectrum produced during inflation in terms of the slow-roU parameters or, alternatively, in terms 
of the scalar spectral ns and and the tensor to adiabatic scalar amplitude ratio r. The bispectrum 
represents a specific non-Gaussian signature of fluctuations generated by quantum oscillations during 
slow-roll inflation. However, our flndings indicate that detecting the non-Gaussianity in the cosmic 
microwave background anisotropies emerging from the second-order calculation will be a challenge 
for the forthcoming satellite experiments. 
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I. INTRODUCTION 



Inflation represents a successful mechanism for the causal generation of primordial cosmological perturbations 
in the early Universe [1]. These fluctuations are then amplifled by the gravitational instability to seed structure 
formation in the Universe, and to produce Cosmic Microwave Background (CMB) anisotropies. Since the primordial 
^\ • cosmological perturbations are tiny, the generation and evolution of fluctuations during inflation has always been 
[ studied within linear theory. On the other hand, there exist physical observables, such as the three-point function of 
■ scalar perturbations, or its Fourier transform, the bispectrum, for which a perturbative treatment up to second order 
I is required, in order to obtain a self-consistent result 
f~| . The importance of the bispectrum comes from the fact that it represents the lowest order statistics able to distinguish 
Oh' non-Gaussian from Gaussian perturbations for which odd-order correlation functions necessarily vanish. An accurate 
Q : calculation of the primordial bispectrum of cosmological perturbations has become an extremely important issue, as a 
number of present and future experiments, such as MAP and Planck, will allow to constrain or detect non-Gaussianity 
of CMB anisotropy data with high precision. 
d • So far, the problem of calculating the bispectrum of perturbations produced during inflation has been addressed 
<^ [ by either looking at the effect of inflaton self-interactions (which necessarily generate non-linearities in its quantum 
fluctuations) in a fixed de Sitter background [2] , or by using the so-called stochastic approach to inflation [3] ^ , where 
back-reaction effects of field fluctuations on the background metric are partially taken into account. An intriguing 
result of the stochastic approach is that the dominant source of non-Gaussianity actually comes from non-linear 
gravitational perturbations, rather than by inflaton self-interactions. 

In this paper we provide for the first time the computation of the scalar perturbations produced during single- 
field slow-roll inflation up to second order in deviations from the homogeneous background. We achieve different 
goals. First, we provide a gauge-invariant definition of the comoving curvature TZ at second order of perturbation 
theory. The importance of the second-order comoving curvature perturbation TZ comes from the fact that it allows to 
compute the three-point correlation function for the primordial scalar perturbations - or its Fourier transform, the 
bispectrum - which represent the lowest order statistics able to distinguish non-Gaussian from Gaussian perturbations. 
Secondly, we show that the second-order comoving curvature perturbation is conserved on super-horizon scales, like its 
first-order counterpart^. Third, we obtain the expression for the gauge-invariant gravitational potential bispectrum 
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during inflation, in terms of slow-roll parameters or, equivalently, of the spectral indices of the scalar and tensor 
power-spectra. 

An accurate calculation of the primordial bispcctrum of cosmological perturbations has become a crucial issue, as a 
number of present and future experiments, such as MAP and Planck, will allow to constrain or detect non-Gaussianity 
of CMB anisotropy data with high precision. 

The plan of the paper is as follows. In Section II we write the perturbations of the metric for a spatially flat 
Robertson- Walker background up to second order and we derive consistently the fluctuations of the energy- momentum 
tensor of a scalar field. In Section III we demonstrate how to find a second-order gauge-invariant definition of the 
comoving curvature perturbation. Section IV will be devoted to the perturbed Einstein equations up to second 
order in the metric and in the inflaton fluctuations. We shall explain how to derive the evolution of the curvature 
perturbation on large scales during a period of cosmological inflation by performing an expansion to lowest order 
in the slow-roll parameters. Finally, in Section V we draw some concluding remarks relating our findings to the 
gauge- invariant gravitational potential bispectrum which is the main physical observable which carries information 
about the primordial non-Gaussianity. We also provide two Appendices where the reader can find various technical 
details. 



II. PERTURBATIONS OF A FLAT ROBERTSON- WALKER UNIVERSE UP TO SECOND ORDER 

We first write down the perturbations on a spatially flat Robertson- Walker background following the formalism of 
Refs. [6,7]. We shall first consider the fluctuations of the metric, and then the fluctuations of the energy- momentum 
tensor of a scalar field. Hereafter greek indices run from to 3, while latin indices label the spatial coordinates from 
1 to 3. If not otherwise specified we will work with conformal time r, and a prime will stand for a derivative with 
respect to r. 



A. The metric tensor 



The components of a perturbed spatially flat Robertson- Walker metric can be written as 

500 = -a' (r) (l + 2<^W+. 



1 



,(2) 



9ij = a^(-r) 



(l_2V;W-V('))(5y 



^(1) 



1 



(2) 



(1) 



The standard splitting of the perturbations into scalar, transverse {i.e divergence-free) vector parts, and transverse 
trace-free tensor parts with respect to the 3-dimensional space with metric Sij , can be performed in the following way: 



-Ar) 



(2) 



xg^ = A,x(^^+a.xr+a,xr^ + x?, (3) 

where (r) = (1), (2) stand for the order of the perturbations, Ui and Xi are transverse vectors = d^Xi'^ = 0)) Xij^ 

is a symmetric transverse and trace-free tensor (^^Xij"* = 0, x^'^P = 0) and Dij = didj — (1/3) 5ij d'^dk is a trace-free 
operator ^. Here and in the following latin indices are raised and lowered using S^^ and Sij, respectively. 
For our purposes the metric in Eq. (1) can be simplified. In fact, first-order vector perturbations are zero in the 



^Notice that our notation is different from that of Refs. [8,9] for the presence of Dij, while it is closer to the one used in 
Refs. [10,11]. As far as the first-order perturbations are concerned, the metric perturbations 'ip and E of Refs. [8,9] are given 
in our notation as ip = f/)*-^^ -|- (1/6) 9i 9*x^^^ and E = x'^V^i respectively. However, no difference appears in the calculations 
when using the generalized longitudinal gauge in Eq. (28). 
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presence of a scalar field; moreover, the tensor part gives a negligible contribution to the bispeetruni. Thus, in the 
following we can neglect and Xij^ ■ However the same is not true for the second order perturbations. In 

the second-order theory the second-order vector and tensor contributions can be generated by the first-order scalar 
perturbations even if they are initially zero [7] . Thus we have to take them into account and we shall use the metric 

goo = -a\r) (l + 2</,W + , 
90i = a'{T) (d^a^W + l-d^J'^ + lcoP 
9ij = a^(r) 



The controvariant metric tensor is obtained by requiring (up to second order) that g^i^g^ = <5/n and it is given by 



(4) 



r = 



a-\r) 



(5) 



Using g^iy and g^'^ one can calculate the connection coefficients and the Einstein tensor components up to second 
order in the metric fluctuations. Their complete expressions are contained in Appendix A. 



B. Energy-momentum tensor of a scalar field 

We shall consider a scalar field <^(t, x') minimally coupled to gravity, whose energy- momentum tensor is given by 

(6) 



where V{(p) is the potential of the scalar field. A successfull period of inflation can be attained when the potential 

V{ip) is flat enough. 

The scalar field can be split into a homogeneous background (Po{t) and a perturbation S(fi{T,x^) as 



ip{T, x') = ipoir) + Sip{T, x') = ifoir) + 6^^^ip{T, x') + ^6^^^<f{T, x') , 



(7) 



where the perturbation has been expanded into a first and a second-order part, respectively. Using the expression (7) 
into Eq. (6) and calculating = g^^T^v up to second order we find 



where T^^°^ corresponds to the background value, and 



(8) 

(9) 
(10) 



2 



a' - 2 
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(11) 

(12) 



(13) 
(14) 



5(2) = 4 



lr(2) 5^ 2 

2 



dip 



(15) 



(16) 



with Vb = y((po)- A comment is in order here. As it can be seen from Eq. (5) and Eqs. (10), (12), (14) and (16) 
the second-order perturbations always contain two different contributions, quantities which are intrinsically of second 
order, and quantities which are given by the product of two first-order perturbations. As a consequence, when 
considering the Einstein equations to second order in Section IV, first-order perturbations behave as a source for 
the intrinsically second-order fluctuations. This is an important issue which was pointed out in different works on 
second-order perturbation theory [12,13,7] and it plays a central role in deriving our main results on the primordial 
non-Gaussianity. 



III. THE SECOND-ORDER GAUGE-INVARIANT COMOVING CURVATURE PERTURBATION 



Let us now focus on the primordial cosmological perturbations produced during a period of inflation driven by 

the inflaton ip. The fluctuations of the inflaton produce an adiabatic density perturbation which is associated with 
a perturbation of the spatial curvature ip. The density/curvature perturbation can be defined in a gauge-invariant 
manner as the curvature perturbation TZ on slices orthogonal to comoving wordlines. At first order, the comoving 
curvature perturbation is given by the gauge-invariant formula [8] 

7^(l) = ^(l) + 2l5(l)<^, (17) 

where Ti = a' /a is the Hubble rate, ■i/'^^-' is the gauge-dependent first-order curvature perturbation and S^^'^ip is the 
first-order inflaton perturbation in that gauge. An important feature of TZ^-^^ is that it remains constant on super- 
Hubble scales while approaching the horizon entry. 

We want to flnd the gauge-invariant comoving curvature perturbation TZ up to second order. Thus we expand the 
curvature perturbation ^ and the inflaton perturbation Sip as tp — 'ip'-^'' + ^ ^^^^ and Sip = S^^^ip + \ S^'^^ip, according 
to Eqs. (1) and (7). An infinitesimal coordinate change up to second order induces a gauge transformation of the 
metric and the scalar field perturbations [6,7]. In particular for a second-order shift of the time coordinate 

---^?i) + ^(«i)-^?2)) (18) 

ip^"^^ and S^'^^ifi will transform as [6,7] 
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= (5(2)^ + ^0^^ (^^Xn) + ^Xn)' + 2'5(^ V) + ^oC['2) • (20) 
We find that the gauge-invariant comoving curvature perturbation TZ = TZ^^^ + ^TZ^"^^ is provided by 

i(vw' + 2Wv(i) + w5W<^'m)' 1 . 



7^ = 7e(i) + 1 



,(5(2) 



+ 



6 



(21) 



We devote the rest of this section to illustrate how to find such a quantity. 

Let us consider the first two terms in Eq. (21) and how they change under the time coordinate shift of Eq. (18). 
According to our perturbative expansion 



V^o 2 

and using Eqs. (19) and (20) this quantity transforms as 



n- 



Xi2) 



^0 



where we have set T = V^^^' + 2WV^^^ + W^^^VVvo- Notice that 



T = T + ^l^ 



n 



(22) 



(23) 



(24) 



since the usual first-order transformations for ip^^^ and S^^^ip are ip^^^i — -ip^^^ — "^C^i) and <5(i)(/? = d^^^ip + VoC^i), 
respectively. Thus Eq. (23) can be written as 



(25) 



Solving Eq. (24) for and using the the first-order transformation a;(i) = w^^^ — ^^'-^^ to express and 9*^°-^^ we 

finally find 



(^() 2n' + 2W-nv'^/^'^ 6 



aia;(i)a'a;(i) , (26) 



which shows that TZ - the combination on the r.h.s. of this equation is indeed gauge- independent. 
Notice that, by replacing ip with the energy density p, we can find in an analogous way a gauge-invariant expression 
for the curvature perturbation on uniform-density hypersurfaces C, which to first order is given by [14] — C^^^ = 
^(1) where Pq is the background energy density. Thus at second order C, = C*^^'' + ^C^^^ is 



n- 



,(5(2 



l(V'(i)' + 2W^(i)+W<5(i)p7p^)' 1 . 



(27) 



IV. EINSTEIN EQUATIONS 



In this section we shall derive the behaviour on large scales of the comoving curvature perturbation at second 
order 7?.(2) introduced in Eq. (21). Our starting point is to calculate the perturbed Einstein equations = 
(^2/2) S^^^'T^l in a generalized longitudinal gauge. Here = SttGn and is the energy-momentum tensor of the 
inflaton field. From the Einstein equations in this gauge we shall pick up an equation for a single unknown function 
- the potential (/i(2) - in a way similar to the procedure used in Ref. [8] to isolate the equation of motion for (f>^^^ in 
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the longitudinal gauge ^. As in Ref. [8], but at the second-order level in the perturbations, an equation linking S'^'^^ip 
and (j)^"^^ holds, so that it is possible to obtain an explicit expression for and hence for the curvature 7?,^^', once 

the equation for (/)*^^) has been solved. Indeed there are many differences with respect to the first-order case, as it will 
be evident to the reader from the details of the calculations which follow. The main difficulties arise due to the fact 
that - contrary to the first-order perturbation theory - also vector and tensor contributions are present, and to the 
fact that the two scalar potential (^^^^ and i/"^^^ differ even in the longitudinal gauge for the presence of source terms 
which are quadratic in the first-order perturbations. 

A. Einstein Equations in the generalized longitudinal gauge 

Up to now we have not choosen any particular gauge. Hereafter we will work in a generalized longitudinal gauge 
defined as 



5oo = -a2(r)(l+2</.«+0(2)), 
ffoi = , 



9ij = a {t) 



(28) 



One can obtain the Einstein equations in this gauge either by using directly the metric tensor in Eq. (28) or by using 

the more general metric of Eq. (4), where no gauge choice has been specified yet, and reduce the equations to the 
longitudinal gauge only at the end. We have performed both the computations to have a cross check for the equations 
obtained. In Appendix A we give the expression for the Einstein tensor in the more general form using the metric 
(4). 

We shall now give the Einstein equations in the longitudinal gauge at first and second order in the perturbations, 
respectively 

<5(i)G°o = K^S^^^T°o implies (29) 

6 ( ^) '^(1) + 6 -i^^^y - 2d, aV^'^ = f ^o'' - V W - -^^'V ^ (30) 
\ a J a \ oip J 

5(i)gO. = K25(i)r° implies (31) 

_ 2 ^a^cA^i) - 2 dii^^^y = -K^ <po' (32) 
5^1) J. = K^J^i^T^ implies (33) 



2^^W' + 4^</,W - 2f-)%W + a,avW + 4-^«' + 2V(^)"-5,5V(^)V\- (34) 
a a \ a / a / 

6(^)G% = ^^^'^ro implies (35) 



3 ^^(^y - + (^)%^'^ + ^<f>^'^ - 12(^)' (V^^))' - 3 {^P^^yy (36) 



''Wc recall that the equation of motion for the potential 0*^^' in the longitudinal gauge is [8] 

^w" _ Q^o^^m + 2 (t^ - ^) 0(1)' + 2(h'- ^n) 4>^^> = o. 



6 



2 gzpr 



= y^^'^n implies 



o 4 \ / 



implies 

2 



(37) 
(38) 

(39) 



f ia,a'=<^(2) ^ ^^(2)' + ^^(2) ^ /^\2^(2) _ 1 a^a/=^(2) + 2^V^')' + V'('^" (40) 
\2 a a \ a J 2 a 

(^(^)') 



(2)" , i(2)" 

x)' +x)' 



<y5 



2 a2y 



In writing the perturbed Einstein equations at second order we have set throughout 0(^) = tl'^^\ since in the longitu- 
dinal gauge at first-order the two scalar potentials are equal in the case of a scalar field, and to obtain Eqs. (36) and 
(40) we have made use of the background relations a" /a = + Ti! and (k2/2) ip'^^ = — Ti'. 
We shall now describe how to isolate the equation for the potential ^(2) _ \Yg ^gg the (0 — 0)-component of Einstein 
equations, the divergence of the {i — 0)-component and the trace of the (z — j)-component, both performed with the 
background metric Sij. Notice; that the divergence and the trace operations make the vector and the tensor modes 
disappear from the equations. Thus, we are left with three equations in the three unknown functions, ^(2)^^(2)^ g^^j 

5(2)^. 

Prom the divergence of the (i — 0)-component of Einstein equations it is possible to recover an expression for 



-<5(2)^ = 



9 I 



where 



(41) 



(42) 
(43) 



and A ^ is the inverse of the Laplacian operator for the three spatial-coordinates. The expression (41) and its 
derivative with respect to the conformal time r can be used in the trace of the {i — j) equation to obtain 

3 3 a \ y \a J \ y a 3 3 

+ fV'W'f + A-^a' + 2-A-ia -At^A-y- 2-«2A-i/3 + K2[i (s^^'^^V 
V / a a [2 V / 



(44) 
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Thus a relation between ip^"^^ and (^^^^ follows from Eq. (44) 



^(2) = ^{2) _ ^-1^^ (45) 

where 7 stands for three times the R.H.S of Eq. (44). Eq. (45) shows that the two scalar potentials ip^"^^ and (j)^'^^ 
differ for quadratic terms in the first-order perturbations, as anticipated above. 

Using Eq. (45) we are now in the position to express the other two unknown functions ^^^^ and S^^'^ip in terms solely 
of Prom Eq. (41) we finally obtain 

Plugging Eqs. (45) and (46) into the (0 — 0) Einstein equation, the equation of motion for (^^^^ is derived 

0(2)" _ + 2(n-^^ <^(2)' + 2(n'- <t>^^^ = (4?) 

_ (n + A-i/3 + K^A-'P' - 7 + - 2^^ A-iy + A'^" 

+ «2 (-i - ia.<5(^)^aM(^)^ - 2 (^(1))%? - i (<5(i)^)' + 2^(i)5(^)^V^) • 

Eq. (47) is our master equation. Before solving it, let us stress two important points. First, no approximation has 
been made up to now. In particular notice that this equation is exact at any order in the expansion in terms of the 
slow-roll parameters. Secondly, the L.H.S. of Eq. (47) is exactly the same as in the equation for 0^^) in the longitudinal 
gauge at first order (see the footnote 4). However, at second order, the key point is that Eq. (47) for (t)^'^'' is not 
homogeneous, but there is a source made up of terms which are quadratic in the first-order perturbations. Notice 
that, as in the first-order calculation, it is not necessary to use the perturbed Klein-Gordon equation to close the 
system of the evolution equations for the fluctuactions. Nevertheless, we also calculated the Klein-Gordon equation 
at second order in the inflaton and metric fluctuactions; this is reported in Appendix B. 



B. The large-scale curvature perturbation TZ^^^ in the slow-roll approximation 

We shall now solve Eq. (47) in order to obtain the expression for the comoving curvature perturbation TZ^^^ defined 
in Eq. (21). First we rewrite Eq. (47) in cosmic time dt = adr, since it is more easy in this way to recognize the 
slow-roll parameters e = —H/H^ and ry = e — (ipo/HipQ) [1]. During a period of inflation e and t] must be <C 1, but 
only at a certain point we will perform an expansion to lowest-order in the slow-roll parameters. Moreover, where 
possible, we shall neglect some terms which give a subdominant contribution on large scales. 

Using cosmic time, Eq. (47) becomes 

_ lA-^d + 2^h(i-^\ A-1/3 + 4^ A-^/3 - 2^ A-I7 + ^-'l 
a a \ Hifo J a ifo 

- (2(^W^y + 8^o^(V'«)'-0(5W^)'-8^oV(i) (JW^)') , 

where we have replaced the term —7 appearing in Eq. (47) with its explicit deflnition given by Eqs. (44) and (45). 

Notice that in the source on the R.H.S. of Eq. (48) there appears always the combination A^^ (k^/J/o — a/aj and its 
derivative with respect to cosmic time. Let us now calculate such a combination. The quantity a deflned in Eq. (42) 
can be rewritten as 
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Using the equation of motion (30) to express > the quantity /?, Eq. (43), turns out to be 



(49) 



a 2 (fio 
+ ^5fc5(iV9'= (did 



Using the equation of motion for ip^^^ in the longitudinal gauge 
which can be derived from Eq. (32) with (f>^^^ = ip^^\ we get 



(50) 



(51) 



_ a ^ K V'o 



a a 2 \ J 

+ -^duS^^^ipd'^ (did'-. 



- 2di d' 



-did'8^'^^ipdkd''5^^^^^\ 



(52) 



and therefore 



-1 ' _ 



(3 a 
a a 



k2 (^0 ^ ^2 



^ +3F 

2 (^0 ^ 



+ ^A-i \dk5^^^^d'' (did'^^^A +did'6^^^ipdkd''ilj'-'^ 



(53) 



The master equation (48) can be simplified if we drop those terms which are next to leading-order in the slow-roll 
parameters 



- 12K2i75(i)^V'^'Vo + Sk^H^ ( S<-^^S - UH^P^'^W + 18H 
+ 4 



(54) 



12^? ^ 1 



9^(5(1^3'' (ai^'^^i)) +did'5^^^ipdkd^'4)''^^ 



2^A-H + A-i7 

■^0 



(55) 



where wc have used Eq. (51) so that -2AH'^ (l + H/H^^ (?A(i))^ -24ffV^(i)?/'(i) = -24ij (t/;(i))^ - 12K2ifd^(i)(^ 

we have used the expression for ((5(1^)' from Eq. (30), and we also explicitly written the combination (53). 
Let us notice that through Eq. (51) we find 



-12//(^«^)-+6(^«)'-.'0(^(iV)\ 



(56) 
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where the last passage is vaUd to lowest-order in the slow-roll parameters. Thus we finally obtain 



H ipo 



+ 4 



2 ipo \ I 



2^A-i7 + A-i7+ ^^^^-^ 



8 fdid'ij^'^^ 



(57) 



Integrating Eq. (57) we find 
<^(2)+i/</,(2) = _i2if 



+ 4 



+ 6 / UW) + 



2 (^0 
12F 

a2(yjo 



-2 



_^A-i7 + A-i7 



-4 J^(d,i^Wd'^W)dt-J- 



dt + 



8 fd,d'ip^^'> 



(58) 



We are now in the position to calculate the second-order comoving curvature perturbation TZ^'^h Prom Eqs. (21), 
(45) and (46) we get 



24a- 



2^A-H + <^(^)-A-S 



(59) 



where the last term is the part of T?.*-^) in Eq. (21) which is quadratic in the first-order perturbations. Note that it 
is of the order of 0{e^,if) and thus we shall neglect it in the following. In Eq. (59) there appears once more the 
combination A~-^ («;^/3/a — a/a) . Thus using Eq. (53) and Eq. (58) we obtain 



ipo ( (5(1 V 



^0 \ ^0 
2 



12- 



2 12i? 
• 2 



dt 



16F 



iQH r 1 



/«2(^0^ 



• 2 ' 9-2 

.^0 K'^V'o 



V-) dt 



+ 



24H'^ f 1 



^0^ 
16H 



y"^A-i [3^5(1^3'= (^iSV^'^) +aia^,5(iV5i5y^] '^^--^^ J ^ 



i6i? r 1 



I6H (po f 1 



-2^^A-7-A-7, 



(60) 
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where the last passage is vahd to lowest order in the slow-roll parameters and we have used the equation of motion (51). 
In Eq. (60) the terms which are not integrated in time can be safely treated in the long-wavelength limit. To lowest 
order in the slow-roll parameters and on large scales, k <C aH, tp^^^ can be considered as constant and 



2 H" " ^0 



^a) = ^L^,(iV = eiI^ (61) 



which can be derived from the equations of motion (51) and (B.4) in the longitudinal gauge for and S^^^ip, 
respectively. On the other hand, from the definition of the comoving curvature perturbation at first order, Eq. (17), 
and using Eq. (61) we can write 7^^^^ = iJ^^^V/'Po to lowest order in the slow-roll parameters and on large scales. 
Thus in these approximations Eq. (61) can be rewritten as 

V;(l)=e7^(^^ (62) 

Performing various integrations by parts in expression (60), using the perturbation equations at first-order and properly 
substracting the contributions in the far ultraviolet, we arrive at the final expression for the comoving curvature 
perturbation 

7e(2) = (,y-3e) (7^(l))%J, (63) 

where 



~lj {i'^^^y dt + {e-r])A-^diR^'^d'R^^K (64) 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have provided a complete analysis of the second-order scalar perturbations during inflation leading 
to the derivation of the gauge-invariant comoving curvature perturbation TZ. 

The comoving curvature perturbation receives a contribution which is quadratic in TZ^^\ The total curvature 
perturbation will then have a non-Gaussian (x^) component. Reminding that the gauge-invariant potential $ (which 
is related to Bardeen's variable [10] by $ = —^h) and the curvature perturbation TZ are related by $ = |7^, the 
following simple relation in configuration space holds 

$(x) = $Gauss(x) + j d^y d^z JC (y, z) $Gauss(x - y)3>Gauss(x - z) + Constant (65) 

which is valid on superhorizon scales and where the constant is such that ($(x)) = 0. Here ^Causs = ^TZ^^^ is a 
Gaussian random field. The non-Gaussianity kernel in momentum space is given by 

X;(ki,k2) = ^(?7-3e) + /x;(ki,k2) , (66) 

D 

where f/c (ki, k2) is directly related to the function I and is first-order in the slow-roll parameters. The gravitational 
potential bispectrum then reads 

($(ki)$(k2)$(k3)) = (27r)3^(3)(ki+k2 + k3) [2 /C(ki,k2)7'$(ki)7'$(k2) + cyclic] , (67) 

where ^^(k) is the power spectrum of the gravitational potential. We can also define an effective "momentum- 
dependent" parameter /nl given by^ 



^Notice that according to our definition of our /nl has the opposite sign of that in Ref. [15]. 
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/nl ~ ^ (ns - 1) + fic (ki , k2) , (68) 

where we have made use of the relation between the spectral index ns = 1 — 6e + 2?7 for the scalar perturbations and 
the slow-roll parameters. In the final expression (68) we have not written possible terms of the order of 0{e^ ,ri'^)AN 
which might be sizeable for certain classes of inflationary models [AN is the number of e-foldings from the time 
at which a given scale crosses the horizon and the end of inflation; for large-scale CMB anisotropics A A'' « 60). 
Whether these terms cancel out, or cquivalcntly whether the curvature perturbation remains constant on super- 
horizon scales also at next-to-leading order in the slow-roll parameters remains an open issue. However, if present, 
terms of order of 0{e^,r]^)AN resemble those found in Refs. [2,16], coming from the self-interactions of the inflaton 
field. The primordial gauge-invariant potential bispectrum leads to a nonzero CMB bispectrum via the Sachs- Wolfe 
effect {AT/T},^ = (1/3)$. 

Deviations from a scale-invariant spectrum can make the primordial non-Gaussianity non- negligible. The possible 
presence of non-Gaussianity in primordial cosmological perturbations is only mildly constrained by existing obser- 
vations [17,18]. Recent analyses of the angular bispectrum from 4- year COBE data [19] yield a weak upper limit, 
|/nl| < 1-5 X 10"^. The analysis of the diagonal angular bispectrum of the Maxima dataset [20] also provides a 
very weak constraint: |/nl1 < 2330. According to Ref. [15], however, the minimum value of |/nl| that will become 
detectable from the analysis of MAP and Planck data, after properly subtracting detector noise and foreground con- 
tamination, is about 20 and 5, respectively. These results imply that detecting non-Gaussianity at the level emerging 
from our second-order calculation will represent a challenge for the forthcoming satellite experiments. 

Note added: While revising this paper, a related work by Maldaccna, astro-ph/0210603, appeared, where the three- 
point correlation function for the curvature perturbation is obtained by computing the cubic term contributions to 
the Lagrangian. Our results, where comparison is possible, agree with those in astro-ph/0210603. 
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APPENDIX A: PERTURBING GRAVITY AT SECOND ORDER 

1. Basic notation 

The number of spatial dimensions is n = 3. Greek indices {a, (3, ^,u, ....) run from to 3, while latin indices 
(a, 6, k, ....m, n...) run from 1 to 3. The total spacetime metric g^i, has signature (— , -|-, -|-). The connection 
coefficients are defined as 

1 ap ( dgpj , dgfip dga^ 

'■01 - 

The Riemann tensor is defined as 



Da -pa -pa i -pa -pA pia pA /' A 0^ 

The Ricci tensor is a contraction of the Riemann tensor 
and in terms of the connection coefficient it is given by 

p p) pa _ p) pa I pa per _ pa per ( A\ 

J^Hu — '-'a'- ^iv '-'IJ.'- va ^ aa'- till av fxa ■ \^-^) 

The Ricci scalar is given by contracting the Ricci tensor 

i? = i?';. (A.5) 
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The Einstein tensor is defined as 



G/j.^ = R/xv — -j^diivR ■ (A.6) 



The Einstein equations are written as G^v = k'^T^^, so that = SttGn, where Gn is the usual Newtonian gravitational 
constant. 

In the following expressions wc have chosen a specific ordering of the terms. In the expressions in which two spatial 
indices appear, such as Eq. (A. 18), we have assembled together the terms proportional to Sij. The intrinsically 
second-order terms precede the source terms which are quadratic in the first-order perturbations. The second-order 
fluctuations have been listed in the following order as ''P^'^K Ll)'''^\ x*'^\ x!^'' find ^ij'' respectively. This 

ordering simplifies the analogy between the first-order and the second-order equations and allows to obtain immediately 
the expressions in a given gauge. 

2. The connection coefficients 

In a spatially flat Robertson- Walker background the connection coefficients are 

^00 = — ; ^oj = — ; ^'ij = — ^ij ; (A-T) 

r^o = ^oi = = 0. (A.8) 

The first-order perturbed connection coefficients corresponding to first-order metric perturbations in Eq. (4) are 

5«r°, = <^«', (A.9) 

^a)rO, = 5,</.W + (A.10) 

s^rio = - + 5^(1)' + av^'^ , (a.ii) 

<5(i)rO, = - 2 ^ <;6(i) S,, - d.d.J'^ - 2 ^ <5,, - ^(^)' S,j - ^ D,,x^') + \ A^X^'^' , (A.12) 

(X (X 0/ A 

5Wrj,^. = -V(^)%,- + ^A.x^^)', (A.13) 
5(i)r;.fc = dji,^^) ji - dki^^'^ 7] + sV^'^ 7,-^ - ^ av^^) jjk + \ djDix^'^ + \ dkD'jx^'^ - \ d'D^kx^'^ . (a.m) 

At second order we get: 

5(2)rOo = ^ - 2 0(1)' + aV^'^ 9fea;(i) + -d'^J'^ dkJ'^ + d'^J'^ dkUJ^^y , (A.15) 
S^^K^ = 4>^'^ + \ 7(5^'^^'^ + ^) - 2 </>('^ di - 2 ^</.(i) - i,^^^ diW^^^ (A.16) 

^(2)r^o = i9^<^(2) + i +a;^(2)) + \( d'J-^)' + + 2^,^^) d' cf>^^^ - cl>^^y d'J^) (A.17) 

2 2a 2\ \ / J 

a a 

5(2)rO = f _ ^^(2) _ 1^(2)' _ ^^(2) ^ 2<^W^W' + 4-0(1)^(1) + d'^J'^dki^^'^ (A.18) 

■'Va 2 a a \ / a 

- ^^'^..(i) a.a;(i)) <5., - ld,dM'> + \ (A,x(^^' + a,xf'' + a.xr'+(x!f)'^ 
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S^'^^k = I {-dj V^'^ S\ - dk V^') 5' J + d' V(') 5jk) + \ {djD\x^'^ + dkD'jX^'^ - dWjkX^'^) (A.20) 

+ I d,d, x^(^) + i {a, xT + du xf - xS) - \ [d'-'"' + 

+ ^(i)Va;(i)5,-, + V^i) {djD\x^^^ +dkD'^x^^^ -d'DjkX^^^) + dj^^^^ D\x^^^ 



3. The Ricci tensor components 

In a spatially flat Robertson- Walker background the components of the Ricci tensor R^^ are given by 

-Roo = — 3 — + 3 ( — ) ; Roi = ; 
a V a / 



Rij — 



a" /a'\2 

a \ a J 

The first-order perturbed Ricci tensor components are 

a a a 

a V « / « 2 



S^'^Rij = 



a a a \0' J " V'*/ 



(A.21) 
(A.22) 

(A.23) 

(A.24) 

(A.25) 



CI CI \ Qj J A 

a 2 ■' 2 2 

At second order we obtain 

} rt } ) rt ) 7/1 } 



2 a 



2'"^ 2a"^ 2^ 2a 2 

6^ 0(1)0(1)' - aVi)5fc0(i) - 30(i)>(i)' + 2^(1)5^5^ 0(1) - a^V^i) a'=0(i) 
a 

+ 6 V,(i) ^(1)' + 6^(1) V^')" + 3 (^^(1)')' - 0(1)'^^ a;(i) + d'^uj^^^ dk<f>^'^ 



(A.26) 
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a V a / a a 

a 

n n 



(1) 



.(2) 



(A.27) 



^ (X ^ ^ Gj ' Qj 

a V a / a 

a 

- V(i)"aiwW - 5^^w'aiwW - isV^i) zj^'^x^D' + ^^'^ dk D^x^^y + ^(^y dk d\x^'^ 

- IdkD'-^x^'^ D^a^^y - \d>^-^x^^^ dkD^a^^y + \D'^"^x^^y d^D^^X^') + \d^^x'-') d^DrnkX^^y . 

The purely spatial part of 6^^^TZfj,i, is very long, and for simplicity has been divided into two parts, a diagonal part 
which is proportional to 5ij, and a non-diagonal part Rff. 



_ ^^(2) + is.SV^') - +4((^Y + ^) U'A' + 4^<^W<^W' 

a 2 2 a \\ a / a J \ y a 

a a \\ a / a J 



(A.28) 



+ 10 



D' 



a 

a V a / a a 

a 

- Z^'^'-x*^^ + i^rx^^^ + -d'u;^'^ d^D^x^^) - ^-0^^^^) D.^x^^y 

Qj (X Zi Qj 



(A.29) 



--a 



,(2)' 



^'^^'^'l + 1 + t) (^^.x^^) +5.xr +5,xf +xgo 

+ l^[D.jx'''' + d,xf + djxr'' + {xS')') + IdAD'^.x'^^ - ia,a'A,x^^) 



-^9,^X^ + UD,,X^^y' + d.xf +d,x 



+ 2,^(1) a,aj(/)(i) - dj(i>^^^ d,i^^^^ - d^cj)^^^ djtp^^^ + 3a,v^^^ajV^^^ + 2^p'^^'> diOj^p'^^'^ 
+ 4^0(1) + (j}(^y didjuj^^'^ + 2^5.(1) Si a^w^i)' + Si av(i) Sia^w^i) 
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a a 

CI d Qi ^ 

+ ^^«'a,x(^^' + ^afcV(^M.i?^_,.x(i) + ^afe^(^M,i^^,x(i) - ^a,v«a'^A,;^« 

+ ^d,dju^'^D%x^'^' - ^afca^a;WAiX«' + ia'^c.«9,AiX«'+ ^a'=a;Wa,A,x^^)' 
-a'=a;(^)afcA,x^i)' + ia'=a;(i)'a,A,x('^ + ^ 5, AiX^^' " ^ a'^u;^^)' ^.AiX^'^ 

+ ^5'=c.(i)a,A,x^^^ + ^d'^J'^ d,D,,x^'^ - ^d'J') d,D,,x^'^ - ^-d,d^J'^ D^,x^') 

a Qj (X Qi 

+ \d^DijX^^^ dkD'^^x^'^ - IdkdiDmjX^'^ D^^x'-'^ - IdkdjDmiX^'^ D^-^d) 
+ IdkdmDijX^'^ D^^x^^^ + ^£''=™X^'^ didjDkmX^^^ + Ja^Z'^'^x^'^ a,-£>„fcX^'^ • 



4. Ricci scalar 

At zeroth order the Ricci scalar R is given by 



R=\-. (A.30) 



The first-order perturbation of R is 



5{^)R = J_ (-Q^did'uj^^^ - 2did'J^y - 2a,av(i) - e^^^)" - 6-<^(i)' - is-^^^)' (A.31) 

\ a a a 

At second order we find 

5(2)i? = -a.aV^2) _ 3^,^(2)' _ 6 ^^(2) ^ 2 9^9^^') - 9-v(')' - sv^')" - a^a^w^^)' (A.32) 

a a a 

+ -dkdiD^'x^^^ + 24— f<^«f + 2afc</)Wa*=</)W + 4</.waiavw 

a 2 a V / 

+ 24-</)W0W' + 6 + 36 ^(1)' + 2afcV'^i)5V^i) - 4V«5ia>W 

a a 

+ 12,^W^W" - 12^WV(1)" - 36-V(')VW + 6afcV'i)5V'i) + 16 V'^') 9^ ^V'^) 

a 

+ 6 -9*=a;W afe(/.« + 12-^Waia^u;« + 4 0(1)9^9^1)' + 2 0W'ai9Vi) 
a a 
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a 

+ 2afca,a;(i)V'=x(^) + 6^9,a,a;(i)i?'=^x(') - D^'x^'^ D^.x^^y - ^ iJ^'^x^D'i),^^')' 

a 4 

_3^2)'fc^(i)2).,^(i)' _ 2dkd'DmiX^'^D'^"^x^'^ + 9^ 5' AmX^'^ i^^X^'^ 
a 



5. The Einstein tensor components 



The Einstein tensor in a spatially flat Robertson- Walker background is given by 



\ a 



G\ = G\ = Q. 

The first-order perturbations of the Einstein tensor components are 



V a / a a 2 

a-^ \ a 2 / 



a o V a / a 

a 2 J ■' 



2 - J- ■ 2 " 2 
The second-order perturbed Einstein tensor components are given by 



a'^ \ \a J a a 4 

- 12 (^^y (,^(i))'-12^<^«V(i)'- SSiV^i^aV^i) - 8^(1) a^aV^i) + 12^V^i)V^i)' 

V / a a 2 a 

2 2 a a 

-2dkUJ^^^ d'^^^^y - 2V.'(i''a,a*a;(i) - (j)^^^ did^ D\x^^^ - 2i)^^^dkd' D\x^^^ 

+ dkdii^^^^ D^'X^^^ - 2-didkU^^^^ D'^x^^^ - 2-dkOJ^^^diD'^x^^^ - dkOJ^^^ d'D\x^^y 



(A.33) 

(A.34) 
(A.35) 

(A.36) 
(A.37) 
(A.38) 



(A.39) 
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8 8a/ 

^(2)6-^0 = If-aV^'^ + + -dkD'^'x^^^' + -did'x'^^y - -didW^^^ - — a*a;(2) (A.40) 

V a 4 4 4 a 

a Va/ \a/ a a 

a V a / 

a 

V a / a a 2 

(X ^ 

- dki^^^y D'^'X^^) + \d'^J^)D\x^^^" + -d'^uj^'^ D\x^'y - Af-YdkLJ^'^D^'^x^'^ 

+ 2^afca;(i)i^''=X« - Id.D'^-^x^'-'D^^x^'y - Id^D^/^y D'^"^x^'^ 

+ ^d'DmkX^'^D'^^X^'y + ^d'DmkX^'y D'^^^x^'^ - In^'^X^'^drnD^lx^^y), 

5i^)G\ = 1 ( - - - D^.x^^^' - \d, d'^y + ^9, d'u^f^ + 8^0(i)a,<^(i) (A.41) 

\ a 4 4 4 a 

+ 40(i)a,^(i)' + 2^0(i)'a,0(i)- 4^(i)'a,^(i) - 4^(i)a,7/>(i)' + y.^^^) a'c^w 

a \ a J a 

+ ^dk^l^^^^ D\x^^y - ^l,^^y dkD\x^^) - dki^^^y D\x^^^ + O^J^^ dkO^ D'^^x^^^ 
-2^d'^u;('^D,,x^'^ + (^)Vu;(i)A.x(^) + d'^u;^'^ d^d, D^lx^'^ 

- ia-a;(^)a,a^A™X^^) + Id^D'-x^'^ D^^x^^y + ^dM^x^^y D''"^x^'^ 

- ld,D^,x^'^ D'^-x^^y - ld,Dr^,x^'y D'^-x^'^) , 

si2)G<^] = + -ct^^^y + 2 - <^(2) _ ( ^) \m _ 1 a. a^^,(2) + ^(2)" (a.42) 

■^0^x2 a a \a J 2 

+ 2^V^2)' + ^a.ai^(2) + la.av(2)' _ \d,d,D^\^^) + 4f^)V0«)' 
a a 2 4 \a/\/ 

_8- - 8-0(i)<?)(i)' - - 20(i)a,5'0(i) - 4 0(1)^(1)" 

a V / a 

-2</,WV«' - 8-(/)W^«' - 2afc^W5V^i) - 4V«ai5V^i) + f^W")' + 8-^(1)^^1)' 
a V / a 

+ 4^(1)^(1)" + 2^«aiav^i) - (/-Wa^a^wW - 2</.waia'a;«' - 2 - ^fc^w a'=<^« 

a 

_ 4-0(1) a^a^c^W + - f-f SfeO^W a'^a;^ + 2-5,a;WaVi)' 

a 2 a V a / a 

- -aiS'^wa^a^wW + - a^s'^wW a^a^^w + 4-^«aia'a;« + 2v(i)aia'a;«' 

22 a 
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a a 
+ ^dkd' DmiX^'^ D''"^X^'^ - ldid'DmiX^'^D"^'x^'^ + ^Omd" Di^x^'^ D^^x^'^ 

8 a y ■' 

s(2)Qnd^ = J_ r _ i d'djcp^^') + 1 d'dji>^^^ - -d'djj^^ - \ d'd^j^y -\- (d'jp + djuj'^^A (A.43) 

■' a I 2 2 a 2 2 a \ ■' ) 

a a 

-2^<A(i)i5^x(i)' - ^0(i)'d^x(i)' - 0(i)i?>(i)" + 1^^0(1) a^iP^.x^i) + ia,0(i)a,i)'=^x« 

- la,0«a'=D^x(i) + a,9,0(i)7^'=^x<^) + ^V^^^^'^^^-x^^^' + ^(^^''^^^-x'^^ + 2^v(^^'i?^x(^^ 

+ 5«i?*:,x<'^ + 2^^«i5^x^i)' + v(i)D>'(i)" + i^feV^i) a,-D'='x^i) - ^a,4/.(i)a^-i5Vx^i) 
+ 2 V'(i) Sfca*£)'=j-x(i) + 2v(i)afeaj£>'=*x^^^ - 2i)(^'i dud^D'-x*^^^ -did'tp^'^^ D'jx^'^^ + S^D^-x^i) 

-d'^u^'^d^D^.x^^y + Id'j^yd^D.jx^'^ + ^d'j^y djD\x^'^ - Id^j^yd^o^.x^^^ 

+ dkdju^^y D'\^^) + -d^J^^d'DkiX^^) + -9'=a;(i)a,D'fex'i) - -d^J^^dkD^^x^^^ 
+ 2^d,d,co^'->D^>^X^'^ - iD'^^^yDk^x^^y - ld^D^,/'^d,D''"^x^'^ 

- ldjD\x^'^ dkD'^^x^'^ + IdmD'jX^'^ dkO'^'^x^'^ - \dkd'D,^jX^'^ D'^^x^'^ 

- ^dkdjD^x^'^ D'^'^X^'^ + IdkdrnD'jX^'^ + ^D'^-^x^') d'djDkmX^'^ 
+ \d^D"^\'^'^d,D^,x^') - duff-D\J'^D^^x''^ - ^D.^x^^yD^^^'^ 

- ^ £>ya;(i)'VV') - dmdkD^.x^'^D'^'x^'^ + Idmd"" D^jX^'^ D'^'x^'^] , 

where S'^^^G'^' j stands for the diagonal part of J^^^GJ, which is proportional to (5*^-, and 6^^^G"'^^j is the non-diagonal 
contribution. 
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APPENDIX B: PERTURBING THE KLEIN-GORDON EQUATION 

In the homogeneous background the Klein-Gordon equation for the scalar field <^ is 

<^o" + 2-^o' = -|^a' (B.l) 
a dip 

The perturbed Klein-Gordon equation at first-order is 

^(1)/ + 2 -,5(1)^' - di a'<5(i)^ - <^(i) Vo' - 3 V-d) Vo' - di a'wd) ^o' = -<5(iV K - 2 ./.(D (B.2) 
a dip-' dip 

At second order we get 

^ ^ ^ 

a 

+ aVi)5fcw(iVo" + 2 -9^1) ^o' + 9fea;(i)aVi)W + 2 9 5fe V' 

a 

+ 2^a'=a;(i)5fc(5(i)(/p + 5^ ^^^V' + 8"^^'^ dkS^'^^ip - d'^co^'^ diD^x^'^ ipo' 

To obtain the Klein-Gordon equation in the longitudinal gauge of Eq. (28) one can simply set x^^'' = X^^^ = 0) and 
(^(1) = ■^(1). Thus at first-order we find 

SW/ + 2 ^5(i)<p' - d, d^S^'^'ip - 4.^(1) Vo' = -<5(^V S - 2 .^(^^ , (B.4) 
a dip'' dip 



while at second order the equation is 



^ £1 ^ 6£ ^ 



(B.5) 



d^V 2 1.(2) 2 l,.m ^2'9''^' 2 



where we have used the background equation (B.l) and the first-order perturbed equation (B.4) to simplify some 
terms. 
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